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Abstract 

Using the stabilized jellium model in two schemes of 'relaxed' and 'rigid', we 
have calculated the dissociation energies and the fission barrier heights for the bi- 
nary fragmentations of singly-ionized and doubly-ionized Ag clusters. In the calcu- 
lations, we have assumed spherical geometries for the clusters. Comparison of the 
fragmentation energies in the two schemes show differences which are significant in 
some cases. This result reveals the advantages of the relaxed SJM over the rigid 
SJM in dynamical processes such as fragmentation. Comparing the relaxed SJM 
results and axperimental data on fragmentation energies, it is possible to predict 
the sizes of the clusters just before their fragmentations. 

1 Introduction 

Metal clusters are one of the building blocks of the nano-structured systems. The stability 
of these systems is of great importance in nano-technology. One of the mechanisms which 
can destroy the stability of metal clusters is charging. These systems have been extensively 
studiedpQ in the context of jellium model (JM). In the JM, the discrete ions are replaced 
by a uniform positive charge background of density n B = 3/47r(rf ) 3 in which rf is the 
bulk value of the Wigner-Seitz (WS) radius of the valence electrons in the metal. A more 
realistic version of the JM, the stabilized jellium model (SJM), which was introduced 



1 



[2] by Perdew et. al. in 1990, has improved some drawbacks of the JM. The SJM 
calculations for the fragmentation processes of metal clusters reported so far were 
based on the assumption that the density parameter r s takes its bulk value rf, which 
we call 'rigid' SJM. However, since the surface effects have large contributions in the 
energetics and sizes of small clusters, a more sophisticated usejl] of the SJM is needed to 
predict the correct energetics of the clusters in the study of the fragmentation processes. 
This method, 'relaxed' SJM, has already been used to predict the equilibrium sizes and 
energies of neutral as well as the charged[5J metal clusters (This method is also called 
self-compression method). In contrast to the JM and the rigid SJM in which the r s 
value is borrowed from the bulk system, in the relaxed SJM, the density parameter of 
the jellium sphere assumes a value such that a cluster with a given number of electrons 
and specific electronic configuration achieves its equilibrium state. Comparing the relaxed 
SJM results with the rigid SJM results on the one hand, and the experiment on the other 
hand, provide information on the possible structural relaxations of the clusters in the 
fragmentation processes. 

In this work, we have studied the binary decay processes of positively charged silver 
clusters Ag^ (Z=l, 2) containing up to 100 atoms, in all possible channels, using the 
relaxed and the rigid SJMs. The possible decay channels for singly-ionized Ag clusters 
are 

AgJ+- Ag^ + AgJ, p = l,2,...,JV-2. (1) 

For doubly charged clusters, the decays can proceed via two different processes. The 
first one is the evaporation process 

Ag^+ - Ag 2 +_ p + Ag°, p = l,2,...,JV-3 > (2) 

in which one of the products is neutral; and the second one is fission into two charged 
products 

Ag^+ - AgJjLp + AgJ + , p = 2, 3, • • • , [N/2]. (3) 

In the evaporation processes, the negativity of the difference between total energies 
before and after fragmentation in a specific channel, 

D z (N,p) = E z (N-p) + E°(p) -E Z (N), (4) 

is sufficient for the decay in that channel to take place. In the above equation, E Z (N) and 
E°(N) are the total energies of Z-ply ionized and neutral A-atom clusters, respectively. 
However, in fission processes a negative value for the difference energy is not a sufficient 
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condition for the decay of the parent cluster. This is because the competition between 
the short-range surface tension and the long-range repulsive Coulomb force may give rise 
to a fission barrier. The height of the fission barriers are calculated using the two-spheres 
approximation [6 . The situation in a fission process is shown in Fig. ^ In figure ^ 
the fission of a Z-ply charged iV-atom cluster into two clusters of respective sizes N%, 
N 2 = N — Ni and respective charges Zi, Z 2 — Z — Z\ is schematically shown. Qf is 
the energy release, B c is the fusion barrier which is the maximum energy of the Coulomb 
interaction of two positively-charged conducting spheres, taking their polarizabilities into 
account. Bf is the fission barrier height which is defined as 

Bf = -Qf + B c . (5) 

The coulomb interaction energy, E c , as a function of their separations, d, for two 
charged metal spheres can be numerically calculated using the classical method of image 
charges jSj. Our calculations show that the maximum of the interaction energy, B c , is 
achieved for a separation do > R\ + i?2- The equality applies for equal cluster radii and 
charges. The electrostatic interaction energies of pairs of equally charged {Z\ = Z 2 = 1) 
clusters as functions of the pair separation are shown in Fig. El for different pair sizes. 
The radii of the clusters in this figure are calculated from R = iV 1//3 rf . In Fig. |HJ the 
values of the Coulomb barrier, B c are plotted for different pair sizes but equal charges 
Z\ = Z 2 = 1 in the rigid SJM. As is seen, when both of the cluster sizes are small, the 
barrier is higher than the case when at least one of them is larger. 

The most favored decay channel in evaporation processes is defined as the channel for 
which the dissociation energy attains its minimum value 

D z (N,p*) = min [D z (N lP )} , (6) 

and the most favored decay channel in fission processes is defined as the channel for which 
the fission-barrier height attains its minimum value 

B f (N,p*) = mm{B f (N,p)}. (7) 

2 Total energies of clusters 

The total energy of an iV-electron Z-ply charged cluster is obtained by solution of the self- 
consistent Kohn-Sham (KS) equations [7j in the density functional 8J theory (DFT) with 
local spin density approximation (LSDA) for the exchange-correlation energy functional. 
The SJM energy for a cluster in the LSDA is given by[3] 
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E S m K, n±, n+] = E m [n h n h n+] + (e M (r 8 ) + w R (r s , rf)) / dr n+(r) 

+(Sv) ws (r s ,rf) J dr 6(r) [n(r) - n+(r)] , 
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where 
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and 
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Here, n = n-p + nj, which satisfies / drn(r) = N — Z, and n+ is the jellium density which 
satisfies J drn + (r) = N. 0(r) takes the value of unity inside the jellium background 
and zero, outside. The first and second terms in the right hand side of Eq. Q are the 
non-interacting kinetic energy and the exchange-correlation energy, and the last term is 
the Coulomb interaction energy of the system. Em is the average Madelung energy. All 
equations throughout this paper are expressed in atomic units (h = e 2 = m = 1, the units 
of length and energy are bohr and hartree, respectively). The quantity (5f)ws is the 
average of the difference potential over the Wigner-Seitz cell and the difference potential, 
Sv, is defined as the difference between the pseudo-potential of a lattice of ions and the 
electrostatic potential of the jellium positive background j2]. 

The SJM ground-state energy [Eq. (JHJ)] for a cluster with iV electrons is a function of 
N, r s , and rf. In the rigid SJM, r s takes the bulk value rf, whereas in the relaxed SJM, 
it takes the equilibrium value, f s (N), which for a cluster in the ground state electronic 
configuration, is obtained by the solution of the equation 



where, the derivative is taken at fixed values of iV and rf. The total energies of the 
cluster in the relaxed SJM and the rigid SJM are given by E SJM = Esj M (N,f s (N),rf) 
and E S j M (N,rf ,rf), respectively. 

3 Results and discussion 

To obtain the relaxed SJM properties, we have performed an extensive self-consistent so- 
lutions of the KS equations along with Eqs. (|8|)-()11|). and have calculated the equilibrium 
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r s and the energy values of Ag^ clusters (Z = 0, 1, 2) for different cluster sizes (N < 100). 
To show the main differences in the equilibrium r s values of these clusters, which are 
appreciable for small clusters, we have plotted in Fig. 01 the corresponding values only up 
to iV = 34. As is obviously seen in the figure, the neutral and singly-ionized clusters are 
self-compressed for all values of N. However, for doubly-ionized clusters, the values cross 
the bulk border (i.e., rf = 3.02) at JV = 7. 

Figure shows the relaxed SJM energies per atom in electron volts for neutral, singly- 
ionized, and doubly-ionized silver clusters with N < 34. For comparison, we have also 
plotted the bulk value (e = —7.89 eV) by a dashed line. Using these values of equilibrium 
energies and sizes, the relaxed SJM fragmentation energies are calculated by the Eqs. 

In Fig. EJ we have plotted the dissociation energies of the most favored evaporation 
channels of the singly ionized clusters. We have shown the most favored value of p by p*. 
The solid small square symbols show the most favored values p* on the right vertical axis 
whereas, the corresponding dissociation energies, D l+ (N, p*), are shown on the left vertical 
axis by large open squares. As is seen in the figure, there exist some maxima and minima. 
The maxima of the D z (N,p*) correspond to the closed-shell Ag^ clusters with iV=3, 9, 
19, 35, 59,. . .. These clusters have high stabilities compared to their neighboring sizes. 
On the other hand, the minima correspond to the sizes which decay into two closed-shell 
clusters (for example, Ag^ ^Ag^+Ag 2 ). A negative value for the dissociation energy 
implies that the cluster is unstable against the spontaneous decay. 

In Fig. 13(a) , we have shown the most favored products Ag®« and the dissociation 
energies D 2+ (N,p*) for the decay of Agjf via evaporation channel. It is seen that, here, 
the mean dissociation energy is higher than that in the evaporation of singly ionized 
clusters. That is, here, the number of clusters stable against the spontaneous evaporation 
is larger than that in the singly ionized case. However, evaporation is not the only decay 
mechanism for multiply charged clusters, and they can also decay via fission processes in 
which both fragments are charged. 

Figure [3(b) shows the barrier heights Bf(N,p*) for the most favored fission channels 
Ag^ — > Ag%_ p , + AgJ. . Here, in the relaxed SJM, the Coulomb barriers are calculated 
using the equilibrium sizes of each Z-ply charged fragments i.e., R Z (N) = N 1 ^ 3 ff(N). 
As is seen, for small clusters, the majority have negative barrier heights. That is, most of 
them are unstable against spontaneous fission. However, as N increases, the number of 
clusters with negative barrier heights decreases and beyond a certain size range, all the 
barrier heights become positive. As before, the fragment sizes p* are shown in the left 
vertical axis. 

As mentioned before, a doubly charged cluster can decay both via evaporation and 
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fission. To estimate at what sizes which decay mechanism is dominant, we have com- 
pared the dissociation energies and the fission barrier heights in Fig. [7(c). At small sizes 
(N < 21), the fission process dominates because, the barrier heights for the fission are 
lower than the dissociation energies for the evaporation. However, the competition be- 
tween the evaporation and fission starts at N—21. This competition continues with some 
fluctuations until the evaporation dominates completely. Ag clusters of sizes ~ 1 nm and 
~ 2 nm contain ~ 40 and ~ 300 atoms, respectively. Knowing the fact that, in doubly 
ionized clusters, the induced evaporation is most favored for N > 21, we conclude that the 
ionized clusters of nano-meter size undergo the fragmentation mostly via the evaporation 
processes. 

The rigid SJM fragmentation energies have been calculated by the self-consistent solu- 
tions of the KS equations for the system with the positive background density parameter 
of rf. In Figs. Efa)-(c), the corresponding fragmentation energies of the relaxed and the 
rigid SJM are compared. As is seen, in some cases the energy difference is as large as 0.5 
eV, indicating the significant effect of relaxation. 

Improvements over our relaxed SJM results can be obtained by assuming that, at the 
instant the fragmentation takes place, the the parent and the products are not at their 
individual equilibrium state but that, the r s values of the products are equal to that 
of the parent at that instant. Thus, in the equations used for the dissociation energies 
and the fission barrier heights this fact must be taken into account. The r s value of the 
parent cluster at the fragmentation instant (i.e., at the instant the two pieces are created) 
is something which is between that of the neutral cluster with the same N and that of 
the relaxed Z-piy charged of the same N. That is, r° s (N) < r s < r 2 +(N) for the fission 
of a doubly charged iV-atom cluster. The reason for this fact lies in the probabilistic 
nature of the energy absorption in the ionization process of the neutral cluster as well as 
the energy absorption in the fragmentation process of the Z-ply charged parent cluster. 
If the fragmentation energy is absorbed as soon as the electrons are detached from the 
neutral cluster, then the parent charged cluster does not have enough time for structural 
relaxation due to its charging. Therefore, the appropriate r s value would be that of the 
neutral cluster. On the other hand, if the absorption of the fragmentation energy takes 
place with enough delay relative to the ionization instant, then the parent charged cluster 
has a chance to arrive at its relaxed r s value. These arguments lead us to conclude 
that the experimental fragmentation energies are highly dependent on the details of the 
experimental setup. This freedom in taking the r s values of the clusters in a dynamical 
process such as fragmentation shows the advantages of the relaxed SJM over the rigid 
one. However, in cases where the relaxed and the rigid r s values are more or less the 
same, one would expect that applying the time-saving rigid SJM would lead to the same 
estimates of the fragmentation energies. 
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4 Summary and conclusions 



In this work, we have calculated the dissociation energies and the fission barrier heights 
of the fragmentations of singly and doubly ionized silver clusters in the two schemes of 
relaxed and rigid SJMs for different cluster sizes, N < 100. In our calculations, we 
have assumed spherical geometries for the clusters. In the rigid SJM, the density of the 
positive charge background is taken to be that of the bulk system, n B . However, in the 
relaxed SJM, that density is obtained in a self-consistent way which corresponds to the 
equilibrium state of the cluster. The fission barrier heights are calculated using the two- 
spheres approximation and for the Coulomb interaction of the two metallic spheres we 
have employed the method of image charges numerically. The differences in the results 
of the two schemes show the advantages of the relaxed SJM over the rigid one in that, 
one is able to take into account the structural relaxations in dynamical processes such as 
fragmentations. 
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Figure 1: Fission barrier in the two-spheres approximation. The parent iV-atom Z-p\j 
charged cluster decays into two clusters of sizes N± and N — Ni, with charges Z\ and Z 2 , 
respectively. 

Figure 2: Coulomb energy in electron volts of two singly charged metallic spheres as a 
function of their separation distance. 

Figure 3: Coulomb barrier heights in electron volts for two singly charged metallic spheres 
with different sizes. 

Figure 4: Equilibrium r s values in atomic units for neutral, singly ionized, and doubly 
ionized silver clusters. The horizontal dashed-line shows the bulk value for silver, 3.02. 

Figure 5: Equilibrium energies in electron volts of neutral, singly ionized, and doubly 
ionized silver clusters. The horizontal dashed-line shows the bulk value for silver, -7.89 
eV. 

Figure 6: Dissociation energies, in electron volts, of the most favored decay channels of 
singly ionized clusters, in relaxed SJM, are shown with respect to the left vertical axis. 
The right vertical axis shows the sizes of the fragments in the most favored channels. 

Figure 7: Relaxed SJM results for doubly ionized silver clusters, a)- the same quantities 
as in FiglHlfor doubly ionized clusters; b)- the fission barrier height in electron volts for 
the most favored fission channel as a function of the cluster size. The right vertical axis 
shows the most favored product sizes of the neutral clusters; c)- comparison of the decay 
energies of via evaporation and fission mechanisms. Competition starts at iV = 21 

Figure 8: Comparison of the most favored channel fragmentation energies, in electron 
volts, for the two schemes of the relaxed and rigid SJMs. a)- dissociation energies of 
singly ionized clusters; b)- dissociation energies of doubly ionized clusters; c)- fission 
barrier heights for the fission of doubly ionized clusters into two singly ionized clusters. 
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